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We consider continuous time random walks (CTRW) and discuss situations pertinent to aging. 
These correspond to the case when the initial state of the system is known not at preparation (at 
t = 0) but at the later instant of time ti > (intermediate-time initial condition) . We derive the 
generalized aging diffusion equation for this case and express it through a single memory kernel. 
The results obtained are applied to the practically relevant case of the equilibrated random walks. 
We moreover discuss some subtleties in the setup of the aging subdiffusion problem and show that 
the behavior of the system depends on what was taken as the intermediate-time initial condition: 
whether it was coordinate of one particle given by measurement or the whole probability distribution. 
The two setups lead to different predictions for the evolution of a system. This fact stresses the 
necessity of a precise definition of aging statistical ensembles. 

PACS numbers: 02.50.-r, 05.40.Fb 



1. INTRODUCTION 

Fractional diffusion equations (FDE) nowadays can be 
considered as an established and common mathematical 
tool with applications widely distributed in natural and 
social sciences, biology, and finance. Usually, one suc- 
cessfully uses them to explain anomalous scaling behav- 
ior of some quantities of interest. However, when finer 
details of anomalous transport processes come to a ques- 
tion, such as aging for example, an investigator, equipped 
with FDE but with the logic of classical diffusion, often 
fails to describe it correctly. An intrinsically asymptotic 
character of FDE is more crucial than in classical dif- 
fusion, since it involves nonlocal time/space operators. 
Therefore, some particular tasks require the knowledge 
of the dynamics underlying the FDE and not only their 
final form. In many cases, such dynamics can be modeled 
by continuous time random walks (CTRW) P, 0, Q . It 
allows to follow all derivation steps in detail and obtain 
useful analytical results. In the present work, we will also 
use CTRW as a basis model. 

The problem of aging is one of the fundamental ques- 
tions relevant in different fields of physics. The word 
"aging" essentially denotes a whatever pronounced time- 
inhomogeneous behavior of a physical system. It applies 
to all situations when the result of a measurement de- 
pends explicitely on the time when this measurement was 
performed (i.e. what is the time elapsed between the in- 
stant when the system was prepared and the instant when 
the measurement started). Aging and memory effects 
are intrinsic features of glasses [4, |^, 0] , colloidal systems 
0) 9 1 granular materials [§, diffusion processes in 
random environment [ll|, etc. Measurement setups and 
physical reasons for such a time-inhomogeneous behavior 
can be quite different, e.g. the physical system can age 
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due to some slow internal processes or due to external im- 
pacts. In what follows we consider the CTRW model to 
demonstrate the nature and appearance of aging effects 
which, in this case, can be completely understood. 

Different aspects of aging in continuous time random 
walks were addressed e.g. in [H, [H, [13, [H, [11]. The 
problem discussed in the present article is related to the 
one of Ref. , where the author considered the aging in 
the subdiffusion process generated by a deterministic dy- 
namical system. We generalize the results of the Ref. 14 1 
and explicitly express the generalized diffusion equation 
describing aging CTRW by using a single memory kernel. 
Further, we consider the situation of random walks with 
truncated power-law waiting time distributions, which 
leads to equilibration. We show that the final form of 
the corresponding diffusion equation depends on whether 
this waiting time distribution possesses the second mo- 
ment or not. In the latter case, the convergence of the 
probability density function (PDF) for the concentration 
of particles to a Gaussian is very slow. 

We moreover stress that the results obtained in the 
Ref. [l^ and generalized in the present work are perti- 
nent to a very specific variant of the aging problem. To 
demonstrate this, we discuss a seemingly similar (but in 
reality very different) approach to aging. Although both 
approaches are intimately related, they lead to different 
predictions for the evolution of the system. We show 
that these approaches also correspond to different ex- 
perimental situations and discuss their applicability and 
limitations. 

The structure of the paper is as follows. In the Sec- 
tion [5] we formulate the aging problem corresponding to 
the setup of Ref. [l^, introduce the generalized master 
equation (GME) approach, and derive the corresponding 
fractional diffusion equation for aging walks. We analyze 
the asymptotic form of the aging diffusion equation and 
express it through a single memory kernel. In the Sec- 
tion [3l we apply these results to the equilibrated random 
walks. We derive the corresponding transport equations 
with distributed order fractional time derivative and dis- 
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cuss particularly interesting behavior of the moments of 
the density distribution. In the Section 21 we discuss an- 
other approach to aging, the relevant experimental sit- 
uation, and the corresponding results. The last section 
(Section [S]) is reserved for conclusions. 



2. MASTER EQUATION FOR AGING 
RANDOM WALKS 

2.1. General considerations 

Consider a system of particles created at = . Af- 
ter the system was created, we let it evolve according 
to its internal dynamical laws during an interval of time 
ti > . The initial distribution of particles at prepa- 
ration, n{x,0) is unknown. At the time ti we labeled 
some particles and thus created a known probability dis- 
tribution n{x,ti) of such marked particles (performed 
a measurement). For example, we could irradiate our 
system with light and create excited states of atoms or 
molecules of the interest inside the region of the order of 
the light beam radius or its penetration depth. Equiva- 
lently, we may consider a single particle that has started 
at time t — from an unknown location and at the time 
ti , was detected at the point xi and labeled. The quan- 
tity of interest is n{x, , the density at a given instant 
of time t2 > ti , which we assume to be experimentally 
accessible. Then the following question can be put: If 
we only know the duration of the previous evolution, ti 
(called aging time), and the position of the particle at 
t ~ ti , what can we say about the future evolution of 
the particle's position? How precise can we predict the 
position of this particle at the time ^2 > , and how 
does this prediction depend on ti and ^2 ? In the case 
of many labeled particles, the same questions apply to 
the profile at ^2 which can be obtained as a convolu- 
tion of the concentration profile at ti and the PDF of 
the single-particle displacements. We note that the den- 
sity of points in a configuration space of some system 
can correspond not only to coordinates of some real par- 
ticles but to whatever other coordinates characterizing 
e.g. temperature or magnetic field. In our explanation 
we, however, confine ourselves to the picture of particles. 

The problem of aging in this setup can be considered as 
the "intermediate-time initial condition" problem: The 
system was created aX t — Q but the initial condition to 
the corresponding transport equation is posed at a later 
instant of time t = ti . For normal, Markovian diffu- 
sion this does not change the overall form of the trans- 
port equation. In non-Markovian cases, especially in the 
ones with long enough memory, it does, as we proceed to 
show. This aging problem is schematically illustrated in 
the Fig. [TJ Assume that the density of labeled particles 
at t = ti at point x = xi is known. Starting from this 
point the particles diffuse, and at ^2=^1+^ acquire 
some distribution (not shown). By n{x,ti,t) we denote 
our theoretical prediction for this distribution. In the 



case when the whole profile n{x,t = ti) is known, the 
concentration profile nA{x,t — 12) at t2 is given by a 
convolution of the former one and n(x, ti,t) . 




FIG. 1: The first aging problem setup. The evolution starts 
at t = from an unknown initial distribution. At t = ti 
the particles with a profile n{x,t = ti) are labeled. At 
t2 = ti + t for the diffusion from a single point we obtain 
the probability density n{x,ti,t) and for the full profile the 
density n^{x,t = 12) ■ The dashed line and the grayscale 
shadow in the {x,t) plane are used to guide the eye for the 
diffusion from a single point. 

For simplicity we start with the CTRW model on a 
one-dimensional lattice; generalizations to the continuous 
case and to higher dimensions are quite evident. Consid- 
er a discrete set of sites marked by an index i . By n.i we 
denote the occupation probability of each site (density of 
particles) . After a certain waiting time at the site i , a 
particle can jump to the two neighboring sites i — l,i + l 
with equal probability 1/2. The waiting time distribu- 
tion at a site is governed by the probability density ^P{t) . 
The properties of this function determine the regime of 

CXD 

diffusion. If the mean waiting time, (r) = J Tip(T)dT , is 



finite, the resulting transport process will be a normal dif- 
fusion. However, even in this case some aging effects can 
still emerge, with the only exception of the exponential 
waiting time distribution [l^ . If the mean waiting time 
is infinite, we are in the situation of anomalous subdiffu- 
sive behavior, where the aging effects are the strongest. 
It is known that in the subdiffusive regime asymptotic 
transport equations have a form of the fractional diffu- 
sion equation, where fractional differentiation appears in 
the temporal part. It is easy to see that in the case of 
the fractional time derivative, the solutions of the corre- 
sponding equation do not posses the semi-group property. 
This is a strong indication that the choice of the starting 
point or an intermediate initial condition may noticeably 
affect the following dynamics. 

To consider this point in detail, we need to derive 
the transport equation which adequately describes ag- 
ing. Our approach is based on the generalized master 
equation and is similar to the phenomenological deriva- 
tion of the diffusion or Fokker-Planck equations using a 
combination of a continuity equation and the equation for 
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currents [18l |. The continuity assumption contains essen- 
tially two balance conditions guaranteeing the probabili- 
ty conservation: a local one (giving the balance between 
the probability gain and loss at one site) and the one for 
transitions between the two sites (representing particle 
conservation during the jumps). A balance equation at 
each site reads 



(1) 



where j~{t) is the loss current, i.e. the probability for 
a particle to leave the site per unit time at time t , and 
j^{t) is the gain current at a site. 

A particle arriving to the site i at time t comes either 
from the left or from the right. Probability conservation 
for transitions between sites then reads 



(2) 



By combining ([T]) and ^ we get a continuity equation 



(3) 



where f{z)dz is the probability to make a jump within 
the time interval between z and z + dz , so that f(z) 
is the time-dependent density of steps. This can be pre- 
sented in the following form: 



/(^) = ^/„(. 



(6) 



ji=0 



Here /„(i) is the probability density that it is exactly 
71-th jump that takes place at time t . This one is given 
by an n -fold convolution of the waiting time probability 
density V(^) with itself. Under Laplace transform with 
respect to z Eq.® reads: 



1 



(7) 



By indexes p, s we will denote the Laplace components 
corresponding to t and ti , and by k the Fourier compo- 
nent, respectively. The Laplace transform of '0'^^ (r, fi) 
with respect to ti is then: 



According to the waiting time distribution, the loss cur- 
rent at time t is connected to the gain current at the 
site at all previous times: the particles which leave the 
site i at the time t (making a step from i to one of its 
neighbors) either were at the same site i from the very 
beginning, or came there at some later time < r < i . 
A probability density to make a step at time t when 
arriving at r is given by the waiting time distribution 
'0(i — r) . Then for the loss current j~(t) we can write: 



Jo 



T)j+{T)dT. (4) 



However, for the aging initial conditions the probability 
to make the first step after the measurement at ti is dif- 
ferent from iP{t) . Here we postulated it to be indepen- 
dent of the spatial position i and denoted by ip^^H^j ^i) • 
We also assumed that the system was created at the time 
t = , and a measurement was performed at t = ti . 
From here on we denote by t the time elapsed from ti 
and denote by n.^^o = ni{ti) the initial condition for the 
system's further evolution. By using Eqs.(I2]lSl) we can 
rewrite dH) as: 

j-(t) = V'(i)(t,ti)n,,o + f ^{t-T) [n,(r) + j-(r)] dr. 
Jo 

(5) 

Now we have to find the expression for the forward wait- 
ing time distribution of the first step '0*^^^ (t, ti) [l^. The 
forward waiting time t is counted starting from the ob- 
servation point. Let us assume that the jump preceding 
ti (numbered j — 1), took place at time tj_i = z . The 
forward waiting time distribution i/;'^^-' (t, ti) is 



^W(T,ti) 



f{z)il!{ti - z + T)dz, 



V's - lo 



'lp{T')dT' 



The double Laplace transform of this function with re- 
spect to both times t and ti has the following form: 



Ips - V'p 



(8) 



With this information at hand, we take the Laplace trans- 
form of (O with respect to t to obtain: 

From the above equation we find the connection between 
loss current and the occupation probability: 



P'ipp i>p\ti) - iip 



1 - Vp 



(9) 



Now we insert © into the Laplace transform of ^ and 
convert it back to the time domain: 



JO 



ini_i(r) -I- ini+i(T) - n,(r) 



1^ 1^ 



dT 
(10) 



where $(t) and $^^''(t, ii) are defined through their 
Laplace transforms: 



i-i)p 



1 - Vp 



(11) 
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2.2. The second memory term 



T2)) . which thus gets the form 



We note now that Eq. (fTO| can be rewritten in the form 
d f 



d , , d f\. 



(12) 



inj_i(T) + ^n,+i(T) -n,{T) 



dr 



where the memory kernel M{t) is the inverse Laplace 
transform of the function Mp = Tpp/(1 — ipp) , which 
is connected with the density of steps f{t) (hav- 
ing the Laplace representation fp = 1/(1 — ipp)) via 
M{t) = f{t) - 5{t) . Thus, M{t) is the density of all 
steps excluding the first one. Evidently, the asymptotic 
properties of the functions M and / are the same. 

Now we concentrate on the temporal asymptotic be- 
havior of the source term. By using Eqs. ([5]) and PT|) . 
for the double Laplace transform of $(^^(t, ii) we can 
write: 



$(1) = 



1 



p — s 
1 

p — s 



1 



1 



1 -v^. 

1 



1 - v^p 
1 



,s(l - ^Jp) 

1 1 



(13) 



s(l - Vp) 



In order to understand the structure of the above expres- 
sion we use one remarkable property of the Laplace trans- 
form. For any function h , such that its Laplace trans- 
form exists, the following property for double Laplace 
transform holds [iTj : 



[h{t + h)]p^s = — 



(14) 



where s and p are Laplace space variables correspond- 
ing to t and ti respectively. By using this fact we can 
rewrite the expression for <i>*^^^ (t, ti) in time domain: 



$(i)(i,ii) = /(t + ii)-/(t) + 5(t). 



(15) 



where f{t) is the density of steps given by the inverse 
Laplace transform of 1/(1 — ^p) . This one, as we have 
already seen, is connected with the first memory kernel 
M{t) via M{t) = f{t) —6{t) , so that it can be rewritten 
in the form <^^^\t,ti) = M{t + h) + S{t + h) - M{t) . 
However, for whatever aged system {ti > ) the delta 
function vanishes, so that the final result 



(16) 



follows. This is a new expression which connects the 
source term with the memory kernel of the generalized 
master equation for aging continuous time random walks 



d , X 



(17) 



d 
dt 



M{t 



1 / ^ 1 
-n,_i(r) -I- -n. 



l(t) -n,(r) 



dr 



[M{t + h)~M{t)] 



1 



^i(^i) + -jTi'i+iiti) - ni{ti) 



In the more general continuous case, when the jumps' 
lengths are no more discrete but rather distributed ac- 
cording to some probability density g{x) , Eq. (|17|) turns 
to an integral equation of the form: 



dt 
d_ 

m 



n{x, t) 



(18) 



M{t 



+ 00 

T) / g{y) [n{x 



y,T) - n(x,T)] dydr 



-OC 



+ 00 



+ [M{t + h)~M{t)] / g{y)[n{x~y,h)-n{x,h)]dy. 



The discrete case, Eq. ([T7j) corresponds to the choice 
g{x) — [S{x — a) + 6{x + a)]/2 with a being the lat- 
tice spacing. The generalized master equations, Eq. (ll7p 
and Eq. (jlSp , contain the standard random walk part rep- 
resented by the first term and an additional term repre- 
senting the memory on the initial conditions. This mem- 
ory term is expressed through the same CTRW memory 
kernel as encountered in the first term on the right hand 
side and vanishes when either ti = , or iP{t) is an 
exponential. 



2.3. Asymptotic form of the aging equation 

Now we would like to obtain the asymptotic form of the 
above transport equation corresponding to large values 
of X and t (or respectively small k and p ) . Assuming 
Tii (t) to change slowly enough as the function of the site 
number i one can change to the continuous description 
introducing the density n{x,t) with x — ia . The dif- 
ference operators in Eg. p?)) can then be considered as a 
discrete approximation to a Laplacian, so that the corre- 
sponding master equations takes the form of generalized 
diffusion equation with the additional memory term 



d_ 

di 



n{x, t) 



Ydt 



M{t - T)An{x,T)dT 



+ [M{t + h) ~ M{t)] y An(a;, 0) (19) 

The same form follows from the more general continu- 
ous form, Eg. p^ . provided the jump length distribu- 
tion g{x) has a finite second moment ( (a;^) < oo ). The 
corresponding proof follows either the standard Kramers- 
Moyal procedure or can easily be obtained as the small- k 
expansion in the Fourier space. 



5 



For exponential waiting time distribution 
tp^t) = exp{—t / tQ) the memory kernel M(t) — l/tp 
and the equation takes form of the ordinary diffusion 
equation 



|n(a:,i) = |-An(x,0 



(20) 



lacking the second memory term. The combination 
K — 0^/2^0 is the usual diffusion coefficient of the pro- 
cess. Another common choice of the waiting time distri- 
bution (determined by both, real practical situations and 
mathematical convenience) is a power law function, e.g. 
^{t) = 7/io(l + t/toy+'> , with < 7 < 1 , where to 
again gives us the characteristic temporal scale of wait- 
ing times. The expansion of its Laplace transform for 
a small p is: tpp ^ I — T{1 — 7)ioP^ + 0{p) , and the 
corresponding memory kernel M{t) is given by 



M{t) 



1 



r(7)i2;r(i-7)' 



(21) 



in which we recognize an expression proportional to the 
integral kernel of the fractional derivative. After sub- 
stituting this expansion into (jl9p and some algebra we 
obtain: 



d_ 

dt 



(22) 



r(i-7) 



An{x, ti). 



where the generalized diffusion coefficient stands for 
— a^/[2r(7)tQ] . We must note that the second mem- 
ory term in Eg. ip^ contains the Laplacian, whereas the 
corresponding term in the result of Ref. does not. The 
probable reason for that is an insufhcient number of terms 
retained in the asymptotic expansion in Ref. . 



3. EQUILIBRATED RANDOM WALKS 

Let us now consider the system which relaxes to a true 
equilibrium. Such a system would correspond to CTRW 
with the waiting time distribution possessing the first 
moment, (r) = TTp{T)dT , which is however so large 
that the intermediate power-law asymptotics of the func- 
tion is still seen. The examples are a theta-truncated 
power-law ipir) — A/{1 + T)^+"0(r — r) or an expo- 
nentially truncated one ipir) — yl/(l -t-r)^"*"" exp(— r/T) 
(A is the normalization constant), or an exponentially 
truncated one-sided Levy-distribution with the Laplace 
transform ipp = exp[-A{p+ 1/A)" + A/X"] . Here A is 
the parameter with the dimension [T]" and of the abso- 
lute value of unity (the typical step time is set to one), 
and A = (Aa/r)^/*-^""' , where (r) is the mean wait- 
ing time. The equilibrated case is interesting because it 
is experimentally relevant: systems are typically creat- 
ed much earlier than measurements are performed and 



have enough time to equilibrate. The time lag between 
the labeling particles at t = ti and final measurement 
at t2 is however small enough to probe nonequilibrium 
dynamics. The general result for the aging problems can 
be easily applied to this concrete example. The equili- 
brated case corresponds to the limit ii — > oo of aging 
CTRW. Analogously to the aging problem, we illustrate 
the equilibrated problem setup in the Fig. ^ The sys- 




FIG. 2: Problem setup for the equilibrated random walks. 
At t = ti we label some of the particles from their uniform 
distribution and create the profile n{x,t = ti) . Later we 
follow only the diffusion of those marked particles and obtain 
at t2 = ti + t the density profile n^(x, t — t2) ■ 

tem prepared at i = evolves during the time ti , with 
ti — > oo (in practice it means that the system was cre- 
ated a long time before the measurement took place and 
develops a uniform distribution of particles). At the time 
t = ti we select some profile of the particles n{x, t = ti) 
and follow only these particles until the time t2 — ti + t . 

We start from the transport equation ((T9| by taking 
its Fourier-Laplace transform: 



pnk^p-nkih) = - 



■fcVp-V^'H'H^i)' (23) 



where we returned to our earlier notation for the second 
memory kernel $i . In the above equation, the wait- 
ing time distribution has a finite first moment, (r) , and 
therefore — 1 ~ (''") P + o{p) . In addition, the limit 
ti ^ oo should be applied. In this case the expression for 
the second memory kernel ^^{t) is simpler than before. 
Its Laplace transform is given by (see ([8]) and (fTT|) ): 



$(1) = 
p 



1 



1 



p{t) 1 - Ipp 



1 



(24) 



We note that in the time domain we now have 
$(i)(t) = l/(r} - f{t) + 5{t) and that f{t) still 
has a physical meaning of the (time-dependent) densi- 
ty of steps. For a pure power-law CTRW this function 
monotonously decays to zero with time, while for equi- 
librated walks it decays to 1/ (r) . The transition be- 
tween the both regimes takes place at at time tc — (r) . 
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However, this does not mean that the second memory 
function vanishes at times which exceed tc : the differ- 
ence 1/ (r) — M{t) stih can behave as a pure power-law 
in some cases. It is important that the second memory 
function is governed by the subleading term in the expan- 
sion of 1/(1 — -ijjp) , whose behavior depends on whether 
the function ^{t) possesses the second moment (we shall 
call this situation "a sharp cut-off") or does not. 

In the case of a sharp cut-off and for a small p , 
tpp = 1 — {t)p + M2P^ + o{p'^) , where M2 is the sec- 
ond moment of waiting times. Now we have: 

1 1 



1-tpp 



{t)p 
1 

{t)p 
1 

{t)p 



+ A/2p2 + o(p2) 

[l + iAh/ {t))p + o{p)]. 
Ah 



This means that for small p 



$1^) = 1 



M2 



o 1 



+ o(l), 



(25) 



(26) 



so that the function $(^'(t) is integrable, with the inte- 
gral being equal to I—M2/ (r) . This in turn means that 
$i(t) decays faster than and its exact decay form 
depends on how many moments does '4>{t) actually pos- 
sess. We note that the existence of the second moment is 
necessary for such a behavior. The overall fractional dif- 
fusion equation for such a process in the long-time limit 
tends to a normal diffusion equation of the type 



dn(x, t) 
dt 



KAn{x,t) + CS{t)Ano, 



(27) 



where K stands for the diffusion coefficient K = j (r) 
and the constant C is given by C = (1 — M2I (r) )(}} , 
i.e. to the ordinary diffusion equation with some cor- 
rection to the initial condition. This correction vanishes 
only for the Markovian random walk process with ex|)0- 
nential waiting time distribution for which M2 = (t) . 

In the case when the second moment does not exist, 
the situation can be vastly different. In order to present 
this type of behavior let us consider a special example 
of the waiting time distribution for which we can convert 
Eq. (j23p in the Fourier-Laplace domain into a distributed- 
order fractional diffusion equation with a source in space 
and time. Let us take ipp in a form 



ipp = exp 



p{t) 



1 + A(p(r))i- 



(28) 



This function is completely monotonic, which, according 
to the Bernstein's theorem, means that it is a Laplace 
transform of a probability density function. To prove 



this it is enough to note that ^pp has a form 



the function hp > , h. 



V\p=0 



with 



and possesses a 



monotonous derivative (see [201). For p — > , Vp be- 
haves as — ^ — p{t) J i-e. the corresponding PDF 
Tpir) has the mean value of (r) . The parameter A is 
a free parameter of the order of unity which governs the 
precise form of the crossover. 

Provided that p is small enough, we can take 
^p « 1 — p (r) /[I -I- A{p (t))^^"] which corresponds to 
1/(1-V'p) = l/(p(r))-hA(p(T))-" . By substituting this 
expression into (j23p we obtain: 



(29) 



pnk,p-no^k = -X [1 + A(p(r))^ "] A;^»T-fe,p 
+ K{T)A{p{T))-"k^no,k- 



In normal space-time domain this equation corresponds 
to 



dn{x, t) 
dt 



= K 



1 + A (r) 

l-a 



l-a T-il-a 



l\n[x, t) 



K 



A(r) 



r(a) 



-t 



-l+a 



A no, 



(30) 



Note the distributed order derivative ( l-|-Ar-^~" o-D* ~" ) 
on the r.h.s. [l^l, which is typical for systems showing 
crossover behavior, and the extremely slowly dissolving 
of the initial condition! 

The difference between the cases of the equilibrated 
random walks corresponding to the waiting time dis- 
tributions with or without second moment (mirrored 
in the difference of the form of the corresponding gen- 
eralized Fokker-Planck equation) is intimately connect- 
ed with the question whether the mean forward wait- 
ing time (the mean waiting time of the first step 
after the beginning of observations) exists. Indeed, 
looking at the limit of the Laplace transform in the 
first argument of the forward waiting time distribu- 
tion for the equilibrated case ( ti 00 ) we see that 
il>''p'^ ^ (1 — ipp)/ (t) p . Therefore, the first moment 

of v(i)(i) , {h) = /;°^(i)(i,oo)idi = ^(1) is 

p=0 

equal to M2/2(r)^ provided M2 , the second moment 
of the waiting time distribution, exists, and diverges oth- 
erwise. 

One of the interesting peculiarities of the equilibrated 
random walks is found in the behavior of its moments. It 
is known that the second moment of equilibrated CTRW 
behaves diffusively at all times. However, the higher mo- 
ments of the density profile are much more exotic when 
the second moment of iP{t) does not exist. From Eq. ((23|l 
we obtain: 



1 l-^p 



21,2 



p p^ (r) 1 — V'p + ijjpu'^k^ /2 
1 l-ipp 



p P^ {t) l — tpp + a^k"^ 



1 a^k 



2l2 



p{t) 

4)2" 



1 a2"fc2» 



p{t) (1-Vp)"- 



(31) 
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Now we easily find the even moments of tlie density by 
the following formula: = (-I)"(d2»/dfc2n)^^|^^^ , 

We thus have: 



Mo,p - l/p 



2iK^{T) 



M' 



2n,p 



(2Ti)!jf"(r)"-^ 
p2(l-Vp)"~i 



Mo{t) - 1 

M2(0 = 2Kt 
Mi{t)=2AK^ {T)U{t) 

(32) 

M2„(t) = (2n)!if" (T)"-V2„(i). 



All odd moments are zero due to the the symmetry of 
the case considered here. The functions f2n{t) are the 
inverse Laplace transforms of f2n,p — l/p^(l — V'p)""^ • 
For the waiting time distribution with infinite second mo- 
ment considered above, its asymptotic Laplace transform 
is as follows: 



f2n,p 



p 



,11+1 



(33) 



In the time domain it corresponds to the following be- 
havior: 



{r) 



— a{n—l) 



f2n{t) 



V[{n - l)a 



2] 



for t > (t) 



(34) 

This means, for example, that the distribution of the 
particles' positions at an intermediate time is consider- 
ably platykurtotic, and only slowly tends to a Gaussian 
in course of time with a typical transition time (r) . We 
note here that such a slow convergence to a Gaussian is 
typical also for other random walks models with truncat- 
ed power-laws, as exemplified by truncated Levy flights 

MM. 



The GME approach used in Sec. [2] relies on a special 
definition of the aging ensemble. The system was created 
at i = and evolves according to the CTRW dynamics 
up to the time ti , where the first measurement of the 
particle's position (or the the particles' positions) takes 
place. No explicite notion about the initial distribution 
of the particles' positions at t = is assumed or used. 
In any case, we follow up only the particles tagged at 
t = ti . The situation considered in the present section 
corresponds to a different ensemble. 

Imagine now that we have two experiments (probably 
with different initial conditions), which are such that the 
exact measured density in one of them is the same as the 
density of labeled particles in the another one. Will the 
further development of the densities be the same? If not, 
what is the difference? As we proceed to show, the two 
situations lead to different results for the higher moments 
of the distribution. It is worth mentioning that the distri- 
bution of the first exit times after the measuring event is 
the only quantity affected by aging and it determines the 
future evolution of the system for both definitions p3| . It 
depends on the initial conditions through the correlation 
between the particle's position and the forward waiting 
time, as discussed in Ref. 12411. 




4. GREEN'S FUNCTION APPROACH AND 
THE PROBLEM OF AGING ENSEMBLES 

In the previous sections, we considered the results for 
a specific setup of an aging problem corresponding es- 
pecially to a displacement of a single particle labeled at 
t = ti in course of its further temporal evolution or of the 
ensemble of such labeled particles. Let us consider anoth- 
er setup for aging problem. Assume that the full profile 
of the density of particles (or the full probability density 
for a single particle) at t = ti is known (for example, 
we started with a thermodynamically large ensemble of 
particles and were able to measure their density at t — ti 
with high enough precision) or is known in principle as 
it is assumed when calculating, say, the correlation func- 
tions [1^. Does it give us additional information which 
might improve the prediction of the particles' positions? 
What would be the following evolution of this profile and 
how will it depend on the aging time? 



FIG. 3: Tlie second aging problem setup. At t — ti the 
full density profile, n[x,t = ii) is measured. By using the 
Green's function it is possible to propagate it back and repro- 
duce the initial distribution at t = Q , n{x, t = 0) . Then it is 
forwarded until the time t2 = ti + 1 to obtain the prediction 
for the future evolution n'^{x,t — t2) . In contrast to the 
first aging problem, diffusion from a single point x\ will de- 
pend on both, the starting point x\ and the aging time ti : 
n{x,x\,t\,i) . The grayscale shadows and dashed line in the 
{x,t) plane show the diffusion from the initial distribution 
and a single point. 



The difference between the two setups, the one of Fig.[T] 
and the one of Fig. [3l on the qualitative level, can be un- 
derstood already now. In the given point xi at the time 
t = ti , there are particles which have arrived different 
times ago and, therefore, will make their next steps also 
at different times [1^ [2y|. In the first aging setup, we 
virtually ignore the dependence of these exit times on Xi 
(the relative position in the full density profile, if the lat- 
ter is known). Therefore, we assume that in all points 
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the distribution of exit times depends only on ti . In the 
second aging problem, we take into consideration the rel- 
ative position of the point xi . It is clear that particles 
which occur to be found at the wing of the density profile 
typically have made much more steps than those found 
close to the center of the profile. Moreover, the particles 
which had to make a lot of jumps probably have arrived 
at their actual position very recently, i.e. start their last 
waiting period not long ago before ti . This implies that 
the distribution of exit times depends on xi which af- 
fects the diffusion of particles out of this point. Therefore 
corresponding profiles in the Fig. [T] and [21 n{x,ti,t) and 
n(x,xi,ti,t) are different. By repeating the same argu- 
mentation for each point in the profile n{x,t = ti) we 
may conclude that and n*^ will be different as well. 
We will come back to this question later in the text. 

Let us now return to the case when the PDF of parti- 
cles' positions at t = ti is known exactly and note that 
the temporal evolution of the PDF from true initial con- 
ditions to the state at time t is described by the linear 
evolution operator S{t) . If the time evolution opera- 
tor has an inverse defined at least on a set of relevant 
PDFs n{x, t) , one can propagate this PDF back in time 
to t — , thus getting n{x, 0) = no{x) and then forward 
in time up to ^2 > ti to obtain n{x,t2) 27]: 



n{x,t2) = S{t2)S ^{ti)n{x,ti 



(35) 



This defines the linear time evolution operator 
S{t2;ti) = S{t2)S-^{ti) which gives us the PDF at the 
time <2 provided the PDF at the time ti is known. 

Let us first find the explicit form of the S{t2;ti) op- 
erator and then discuss the difference between this ag- 
ing problem, and the first one, considered above. The 
time evolution operator, S{t2;ti) , follows from the ex- 
plicit form of the solution of the transport equation via 
the Green's functions method: it is the integral operator 
containing the Green's functions of the standard CTRW 
equation, G(x, x' ,t) . Let us recall the form of the trans- 
port equation for CTRW without aging, i.e. all particles 
were introduced at t = without history (the first part 
of Eq.dlOl) or Eq. (fT8l) ): 

t 

h{x,t) ^ / ^(t-r) / g{y)[n{x - y,T) - n{x,T)]dydT. 



(36) 

The Green's function is defined as the solution of the 
above equation with n{x,t = 0) = S{x — x') . Then, for 
arbitrary initial distribution no{x) we can write: 



n(x, t) = S{t)n{x, 0) 



G{x,x\t)nQ{x')dx' . (37) 



For the homogeneous situation considered here, the 
Green's function depends only on the difference of its spa- 
tial variables: G{x,x',t) = G{x — x\t) . Therefore, the 
solution can be found as a convolution of the initial con- 
dition with the Green's function. By denoting the con- 
volution operation with *, we can write S{t) = G{x, t)* , 



n{x, t) — S{t)no{x) = G{x, t) * no(x). 



(38) 



The PDFs of the particles' distributions at time ti and 
t2 = ti + t are thus given by 

n{x,ti+t) — G{x,ti + t) * no{x), 

n{x,ti) = G{x,ti) * no{x) (39) 

Under the Fourier transforms, convolutions are changed 
into simple products of Fourier components: 



nkih+t) = Gkih +t2)no^k, 
nk{ti) = Gk{ti)no,k- 

From these two we formally find: 

nk[ti +t) = ■ nk[ti) 



(40) 



(41) 



and now write the time evolution operator 
S{t2,ti) — S{t2)S^^{ti) as an integral convolution 
operator T{x,t2,ti)* with the integral kernel being the 
inverse Fourier transform of Gk{ti + t)/Gk{ti) . 

For the CTRW model its Green's function can be eas- 
ily found analytically |2|], so that also the corresponding 
operator S(t2,ti) can be immediately calculated. The 
Green's function for a CTRW equation (the solution of 
([36l) ) in the Fourier-Laplace space is 



Gk,p — 



p{l- tppQk) 



(42) 



We can now compare the resulting density distributions 
obtained in the two aging problems, namely Eq. (|18p and 
(|¥l]l . To do this we analyze the difference of the two 
densities 6n{x,t) = nP — , with n'^ being the PDF 
predicted by the Green's function approach (j4T]l . and 
being the one given by the by the probabilistic approach 
(ITS)) of Section 12.21 In particular we concentrate on the 
behavior of moments of the density difference defined as 

As an example we consider the following situation. 
Imagine that the system evolved starting from a sharp 
(delta-function) initial condition n{x,t = 0) = 6{x) . In 
this case the particles' distribution at t ~ ti is given by 
the inverse Fourier transform of nk{ti) = Gk{ti) ■ This 
distribution is now used as the intermediate-time initial 
condition at t ~ ti in both approaches, in the one of Sec. 
[2] and in the one of the present Section. Then the results 
of these two approaches are compared. Eq. ([iT|) reduces 
in this case to a simple formula n'j^{t + ti) = Gk{t + ti) . 

Omitting the details of rather tedious calculations, we 
can show that the second moment of the density differ- 
, is equal to zero. However, already the next 



/Sn 



deviates from zero. The final ex- 



even moment, 

pression for the forth moment of the density difference 
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reads: 

(,4x ^ 24 [ (43) 

\ /Sn 7r7r(7)2 J (t-r)i-'r ^ ' 



tl 

J \T J Tl+T 


In Fig.l we plot (a;^)^„ given by (|43|) as a function of 
time for three different aging times, for the case 7 = 1/2 
(the exponent in the power tail of the waiting time distri- 
bution) . These findings show that the approach based on 




t 



FIG. 4: The forth moment of 5n{x, t) — -nP — , as a func- 
tion of time, t, for different aging times ti and 7 = 1/2 . 
Its deviation from zero shows the difference between the pre- 
diction given by Ea.p8[l as compared to the Green's function 
result (|4T1l . 

the backward- forward propagation delivers a result which 
is different from the one given by the approach of Sec. 
[51 It is also necessary to stress, that the aged propaga- 
tor n{x,t,ti) given by (jl8p is not the Green's function 
of aged anomalous diffusion, since it does not reproduce 
the PDF of the particles exactly. The Gaussian situation 
showing no aging is the only one when the both approach- 
es are equivalent (and give the same result): in this case 
Gk{t) oc exp{-Dk'^t) so that Gk{ti+t)/Gk{ti) = Gk{t) . 

Let us discuss some implications of the two approach- 
es. The Green's function approach utilizes more informa- 
tion about the intermediate stage of the system (the full 
knowledge of the whole distribution) and corresponds to 
a different nonequilibrium statistical ensemble than the 
one of Sec. [D On the other hand, the approach of Sec. 
assumes no knowledge about the initial condition, i.e. es- 
sentially starts from the assumption that at t = the 
distribution of particles in space was uniform (which is 
the most reasonable assumption about the state of the 
statistical system provided no additional information is 
given and corresponds exactly to the Jaynes' information 
approach to statistical mechanics |29|). It is clear that 
the homogeneous distribution does not evolve in course 



of the time and stays uniform up to the time t = ti . 
At t — tl we choose some profile n{x,t = ti) from this 
distribution (see the Fig. [5]) and follow its evolution. 

As we have mentioned above, the physical reason for 
being able to predict exactly the dynamics of particles af- 
ter a measurement, is the precise knowledge of their exit 
times. The Green's function approach uses this informa- 
tion only indirectly, just propagating the profile back and 
then forward and automatically provides the correspond- 
ing distributions. This approach is not literally applica- 
ble in a whatever case when the intermediate time density 
distribution is not a full profile resulting from the previ- 
ous evolution of some initial density. There exist, how- 
ever methods which work with microscopic distributions 
directly and can be applied to such cases as well. The 
first one is based on the generalized transport equation 
with microscopic details taken into account [iTj , another 
one considers two-point probability distributions [2^, ■ 
Both of these more general approaches (not restricted to 
the aging problems) give answers corresponding to the 
Green's function approach when starting from a concen- 
trated initial condition (and therefore for a whatever ex- 
actly known one). Moreover, these approaches allow for 
the full solution of the aging problem in the setup of the 
Sec. m where the exit times after the measurement tend 
to be independent of coordinate. The results for the equi- 
librated walks provided by methods of Refs. [13, 113] and 
by that of the Sec.[l]coincide since the information about 
the initial state of the system is forgotten. 

It is important to state that the two ensembles corre- 
sponding to the two approaches discussed in this work are 
not the only two alternatives, but two of a quite broad 
spectrum of possibilities. They provide however the two 
limiting situations: the full knowledge and no knowledge 
about this distribution. There are various situations pos- 
sible, when the distribution of the particles' positions at 
the intermediate time is known to some extent or with 
some uncertainty, in which case the results for aging will 
depend on what exactly was measured at t — ti and 
what the precision of this measurement was. 



5. CONCLUSIONS 

Summarizing our findings, we have derived a general 
equation for the PDF of the particle's position in aging 
CTRW, corresponding to the "intermediate-time initial 
condition" setup, in which the coordinate of the parti- 
cle is measured at some time ti after the preparation 
of the system considered to take place at t = . We 
were able to express the second memory function, de- 
scribing the dissolving of the initial condition, through 
the memory kernel of the CTRW equation. We discussed 
the ensuing forms of the transport equation for the case 
of CTRW with power-law waiting time distributions lack- 
ing the first moment, as well for the equilibrating situ- 
ations, where the first moment exists. We have shown 
that the exact form of the corresponding equation for 



10 



equilibrating walks depends on whether the second mo- 
ment of the waiting time distribution exists or not. In 
the first case, the normal diffusion equation appears, in 
the second case, the evolution is described by the diffu- 
sion equation with distributed-order temporal derivative 
and with a very slowly dissolving second memory func- 
tion following the power law. The asymptotic behavior 
of the higher moments of the density profile for this case 
was also calculated and indicated a slow convergence to 
the Gaussian behavior. 

Moreover, as we tried to show, the problem of aging 
is a very delicate task even in the framework of exact- 
ly solvable model of CTRW. To illustrate this we have 
considered and solved another possible aging setup cor- 
responding to the full knowledge of the PDF of the par- 
ticles' positions at time ti after the preparation of the 



system at t = . The two setups give different predic- 
tions for the evolution of the particles density and corre- 
spond to different experimental realizations. This shows 
that there is still a room for further developments in the 
area. One of those is the question of the prediction of 
the evolution based on an incomplete intermediate-time 
initial condition. This work is currently in progress. 
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